/*
Article Name: Spraying Conflict: Aerial Drug Eradication and Armed Violence in Colombia
Authors: Campos, Nieto-Matiz and Schenoni
Dofile name: Estimations and Descriptive Statistics
Date: 03/2025
*/


use "\Consolidated Yearly Dataset.dta", replace

////////////////////////////////////////////////////////////////////////////////

//# 1. Aerial spraying effect on armed conflict violence (Yearly Estimations)
				
	//# 1.1. Defining variables globals and panel dataset
		/* Time period */ 
		global time_period "year>=2000 & year<=2015"
		keep if $time_period
		
		/* Coca-producing municipality */ 
		global coca_municipalities "coca_hist_2000_2015" // Municipality type (coca or not)
		
		/* Dependent variables */ 
		global dependent_violence_variables "percap_total_events percap_insur_events percap_param_crimorg_events"
		global dependent_spatial_variables "mpc_total_events mpc_insur_events mpc_param_crimorg_events " // Spatial

		/* Independent variables */ 
		global indep_eradication_variables "errad_aerea"
	
		/* Control variables */
		global controls_light "cov_light_mean" // Light control
		global controls_weather "cov_precip" // Precipitation control
		global controls_economic_performance "cov_percap_y_total cov_desemp_fisc" // Economic controls
		global controls_distance_airports "cov_dist"	

		/* Defining panel dataset */		
		xtset codmpio year, yearly
		
	
	//# 1.2. Yearly descriptive statistics - 2000 to 2015
		// cd "$results\Raw tables"
		
		/* Descriptive statistics (mean, sd, min, and max) */
		sum $dependent_violence_variables d_errad_aerea_p90 $indep_eradication_variables errad_manual coca_hist_2000_2015 $controls_light $controls_weather $controls_economic_performance $controls_distance_airports gamazonia gandina gcaribe gorinoquia gpacifica  if year>1999 & year<2016
		
		/* Descriptive statistics (p25, median, p75, p90) */
		tabstat $dependent_violence_variables d_errad_aerea_p90 $indep_eradication_variables errad_manual coca_hist_2000_2015 $controls_light $controls_weather $controls_economic_performance $controls_distance_airports gamazonia gandina gcaribe gorinoquia gpacifica if year>1999 & year<2016 & coca_hist_2000_2015==1, stat(p25 p50 p75 p90)
		
		/* Eradication intensity graph */
		preserve
		collapse (sum) errad_aerea errad_total, by(year)
		line errad_aerea errad_total year if year>1999 & year<2016, xlabel(2000(1)2015, labsize(vsmall) nogrid) ylabel(0(25000)225000, labsize(vsmall) nogrid format(%9.0fc)) mcolor(dknavy) graphregion(color(white)) plotregion(lstyle(foreground)) legend(rows(1) order(1 "Aerial eradication" 2 "Total eradication") size(small) pos(6)) lcolor(dknavy maroon) lpattern(solid dash) lwidth(medthick medthick) /*title("Aerial Spraying Evolution in Colombia", color(black))*/ xtitle("Year", size(small)) ytitle("Eradicated coca crops hectares", size(small))
		// graph export "$results\Figures\Eradication evolution.png", replace
		restore
		
		
		/* Matching Analysis */ // Average variables through period
		preserve
			
			
			collapse (mean) percap_total_events percap_insur_events percap_crimorg_events cov_precip cov_percap_y_total cov_desemp_fisc cov_indrural cov_dist MPI_2005 total_expenditure ln_pop_tot (max) errad_aerea_dummy area_total elevation coca, by(codmpio)
			
			* 1. Estimating propensity scores
				logit errad_aerea_dummy cov_precip cov_percap_y_total cov_desemp_fisc cov_indrural cov_dist MPI_2005 total_expenditure ln_pop_tot area_total elevation // coca omitted due to omitted observations (collinearity)
				predict pscore

			* 2. Estimating the ATET on all our baseline dependent variables
				teffects psmatch (percap_total_events) (errad_aerea_dummy cov_precip cov_percap_y_total cov_desemp_fisc cov_indrural cov_dist), ate nn(2)
				// outreg2 using "$results\Raw tables\M1 Matching.xls", excel
				teffects psmatch (percap_insur_events) (errad_aerea_dummy cov_precip cov_percap_y_total cov_desemp_fisc cov_indrural cov_dist), ate nn(2)
				// outreg2 using "$results\Raw tables\M1 Matching.xls", excel
				teffects psmatch (percap_crimorg_events) (errad_aerea_dummy cov_precip cov_percap_y_total cov_desemp_fisc cov_indrural cov_dist), ate nn(2)
				// outreg2 using "$results\Raw tables\M1 Matching.xls", excel
			
		restore
	
	
	
	//# 1.3. Estimations: Effect of aerial spraying on violence (2000-2015)
	
	* Redefining dependent violence variables names
	global dependent_violence_variables "total_events insur_events param_crimorg_events"

/////////* 1.3.1. Baseline Model */
{
		foreach j of global indep_eradication_variables {
			foreach k of global dependent_violence_variables {
		
				xtreg percap_`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $time_period, fe // Baseline
				// outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes, Control - Light, No)
		
				xtreg percap_`k' `j' Lpercap_`k' $controls_light $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $time_period, fe // Nightime luminosity control inclusion
				// outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lpercap_`k' $controls_light $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Year FE, Yes, Lagged Y, Yes, Controls, Yes, Control - Light, Yes)
		
				xtreg percap_`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $coca_municipalities & $time_period, fe // Coca municipalities only
				// outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, Coca Only, Municipal FE, Yes, Year FE, Yes, Lagged Y, Yes, Controls, Yes, Control - Light, No)
				
				xtreg percap_`k' `j' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $time_period, fe // Excluding lagged outcome variable
				// outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, No, Year FE, Yes, Controls, Yes, Control - Light, No)
			} 
		}
	
			// Coef plots and marginsplots
			xtreg percap_total_events errad_aerea Lpercap_total_events  $controls_weather $controls_economic_performance $controls_distance_airports mpc_total_events i.year if $time_period, fe
			estimates store A
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fA, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(0(20)80)
			// graph export "$results\Figures\1. Marginsplot M1 baseline (overall).png", replace
			
			xtreg percap_insur_events errad_aerea Lpercap_insur_events $controls_weather $controls_economic_performance $controls_distance_airports mpc_insur_events i.year if $time_period, fe
			estimates store B
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fB, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(0(20)80)
			// graph export "$results\Figures\1. Marginsplot M1 baseline (insurgents).png", replace
		
			xtreg percap_param_crimorg_events errad_aerea Lpercap_param_crimorg_events $controls_weather $controls_economic_performance $controls_distance_airports mpc_param_crimorg_events i.year if $time_period, fe
			estimates store C
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fC, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(0(20)80)
			// graph export "$results\Figures\1. Marginsplot M1 baseline (param).png", replace
			
			// combomarginsplot fA fB fC, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))
			// graph export "$results\Figures\Marginsplot M1 baseline.png", replace
			
			coefplot (A, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (B, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(C, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(0(0.001)0.002, glpattern(solid) glwidth(vthin) glcolor(gs15))
			// graph export "$results\Figures\Coefplot M1 baseline.png", replace
}


		****** WITH MANUAL ERADICATION
{
		foreach j of varlist errad_manual {
			foreach k of global dependent_violence_variables {
		
				xtreg percap_`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $time_period, fe // Baseline
				// outreg2 using "$results\Raw tables\M1 `k' Manual Erad.xls", keep(`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes, Control - Light, No)
		
			} 
		}
	
			// Coef plot
			xtreg percap_total_events errad_manual Lpercap_total_events  $controls_weather $controls_economic_performance $controls_distance_airports mpc_total_events i.year if $time_period, fe
			estimates store Am
			margins, at(errad_manual=(0 30000)) atmeans vsquish saving(fAm, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(-20(20)60)
			// graph export "$results\Figures\2. Marginsplot M1 w Manual Eradication (overall).png", replace

			
			xtreg percap_insur_events errad_manual Lpercap_insur_events $controls_weather $controls_economic_performance $controls_distance_airports mpc_insur_events i.year if $time_period, fe
			estimates store Bm
			margins, at(errad_manual=(0 30000)) atmeans vsquish saving(fBm, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(-20(20)60)
			// graph export "$results\Figures\2. Marginsplot M1 w Manual Eradication (insurgents).png", replace

		
			xtreg percap_param_crimorg_events errad_manual Lpercap_param_crimorg_events $controls_weather $controls_economic_performance $controls_distance_airports mpc_param_crimorg_events i.year if $time_period, fe
			estimates store Cm
			margins, at(errad_manual=(0 30000)) atmeans vsquish saving(fCm, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(-20(20)60)
			// graph export "$results\Figures\2. Marginsplot M1 w Manual Eradication (param).png", replace

			// combomarginsplot fAm fBm fCm, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))
			// graph export "$results\Figures\Marginsplot M1 w Manual Eradication.png", replace
			
			coefplot (Am, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (Bm, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(Cm, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_manual) legend(position(6)) coeflabels(errad_manual="Manual Eradication") vertical yline(0, lcolor(gs6)) ylab(, glpattern(solid) glwidth(vthin) glcolor(gs15))
			// graph export "$results\Figures\Coefplot M1 w Manual Eradication.png", replace
		
}
		

/////////* 1.3.2. High intensity dummy */
{
		foreach j of varlist d_errad_aerea_p90{
			foreach k of global dependent_violence_variables {
		
				xtreg percap_`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $time_period, fe
				// outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Year FE, Yes, Lagged Y, Yes, Controls, Yes)
			}
		}
		
		// Coef plot
			xtreg percap_total_events d_errad_aerea_p90 Lpercap_total_events  $controls_weather $controls_economic_performance $controls_distance_airports mpc_total_events i.year if $time_period, fe
			estimates store A1
// 			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fA1, replace)
			
			xtreg percap_insur_events d_errad_aerea_p90 Lpercap_insur_events $controls_weather $controls_economic_performance $controls_distance_airports mpc_insur_events i.year if $time_period, fe
			estimates store B1
// 			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fB1, replace)
		
			xtreg percap_param_crimorg_events d_errad_aerea_p90 Lpercap_param_crimorg_events $controls_weather $controls_economic_performance $controls_distance_airports mpc_param_crimorg_events i.year if $time_period, fe
			estimates store C1
// 			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fC1, replace)
			
			coefplot (A1, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (B1, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(C1, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(d_errad_aerea_p90) legend(position(6)) coeflabels(d_errad_aerea_p90="High intensity aerial eradication (p90)") vertical yline(0, lcolor(gs6)) ylab(-10(10)20, glpattern(solid) glwidth(vthin) glcolor(gs15))
		//	graph export "$results\Figures\Coefplot M1 high intensity.png", replace
}
	
	
	
/////////* 1.3.3. Spatial displacement of violence */
{
		foreach j of global indep_eradication_variables {
			foreach k of global dependent_spatial_variables {

				xtreg `k' `j' L`k' $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
				// outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' L`k' $controls_weather $controls_economic_performance $controls_distance_airports) excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes)
			} 
		}

		// Coef plot
			xtreg mpc_total_events errad_aerea Lmpc_total_events $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store A2
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fA2, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(-10(10)40)
			// graph export "$results\Figures\3. Marginsplot M1 neighbors (overall).png", replace
			
			xtreg mpc_insur_events errad_aerea Lmpc_insur_events $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store B2
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fB2, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(-10(10)40)
			// graph export "$results\Figures\3. Marginsplot M1 neighbors (insurgents).png", replace

		
			xtreg mpc_param_crimorg_events errad_aerea Lmpc_param_crimorg_events $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store C2
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fC2, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(-10(10)40)
			// graph export "$results\Figures\3. Marginsplot M1 neighbors (param).png", replace
			
			// combomarginsplot fA2 fB2 fC2, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the neighboring municipalities' mean number of violent events" "per hundred thousand population.", size(small))
			// graph export "$results\Figures\Marginsplot M1 neighbors.png", replace

			coefplot (A2, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (B2, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(C2, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(, glpattern(solid) glwidth(vthin) glcolor(gs15))
			// graph export "$results\Figures\Coefplot M1 neighbors.png", replace
		
}


		****** WITH MANUAL ERADICATION
{
		foreach j of varlist errad_manual {
			foreach k of global dependent_spatial_variables {

				xtreg `k' `j' L`k' $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			//	outreg2 using "$results\Raw tables\M1 `k' manual erad.xls", keep(`k' `j' L`k' $controls_weather $controls_economic_performance $controls_distance_airports) excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes)
			} 
		}

		// Coef plot
			xtreg mpc_total_events errad_manual Lmpc_total_events $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store A2m
			margins, at(errad_manual=(0 30000)) atmeans vsquish saving(fA2m, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(0(10)50)
			//graph export "$results\Figures\4. Marginsplot M1 neighbors w Manual Eradication (overall).png", replace

			
			xtreg mpc_insur_events errad_manual Lmpc_insur_events $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store B2m
			margins, at(errad_manual=(0 30000)) atmeans vsquish saving(fB2m, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(0(10)50)
			//graph export "$results\Figures\4. Marginsplot M1 neighbors w Manual Eradication (insurgents).png", replace

		
			xtreg mpc_param_crimorg_events errad_manual Lmpc_param_crimorg_events $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store C2m
			margins, at(errad_manual=(0 30000)) atmeans vsquish saving(fC2m, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(0(10)50)
		//	graph export "$results\Figures\4. Marginsplot M1 neighbors w Manual Eradication (param).png", replace


			
		//	combomarginsplot fA2m fB2m fC2m, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via manual eradication") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))
		//	graph export "$results\Figures\Marginsplot M1 neighbors w Manual Eradication.png", replace

			
			coefplot (A2m, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (B2m, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(C2m, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_manual) legend(position(6)) coeflabels(errad_manual="Manual Eradication") vertical yline(0, lcolor(gs6)) ylab(, glpattern(solid) glwidth(vthin) glcolor(gs15))
		//	graph export "$results\Figures\Coefplot M1 neighbors w Manual Eradication.png", replace
		
}
	
	
	
/////////* 1.3.4. Mechanisms exploration */
{
		* Redefining dependent violence variables names
		global mulhom_mechanism_variables "total_mulhom insur_mulhom param_crimorg_mulhom"

		foreach j of global indep_eradication_variables {
			foreach k of global mulhom_mechanism_variables {

				xtreg percap_`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k' i.year if $time_period, fe
			//	outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lpercap_`k' $controls_weather $controls_economic_performance $controls_distance_airports mpc_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes)
			} 
		}
		
		global AoW_mechanisms_variables "total_statev insur_statev param_crimorg_statev turf_violence "

		foreach j of global indep_eradication_variables {
			foreach k of global AoW_mechanisms_variables {

				xtreg `k' `j' L`k' $controls_weather $controls_economic_performance $controls_distance_airports m_`k' i.year if $time_period, fe
			//	outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' L`k' $controls_weather $controls_economic_performance $controls_distance_airports m_`k') excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes)
			} 
		}
		
		
		* Violence against the state
		{
			// Coef plot
			xtreg total_statev errad_aerea Ltotal_statev $controls_weather $controls_economic_performance $controls_distance_airports m_total_statev i.year if $time_period, fe
			estimates store A4
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fA4, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(0(0.5)1.5)
			//graph export "$results\Figures\5. Marginsplot M1 Violence against state (overall).png", replace
			
			xtreg insur_statev errad_aerea Linsur_statev $controls_weather $controls_economic_performance $controls_distance_airports m_insur_statev i.year if $time_period, fe
			estimates store B4
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fB4, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(0(0.5)1.5)
			//graph export "$results\Figures\5. Marginsplot M1 Violence against state (insurgents).png", replace

		
			xtreg param_crimorg_statev errad_aerea Lparam_crimorg_statev $controls_weather $controls_economic_performance $controls_distance_airports m_param_crimorg_statev i.year if $time_period, fe
			estimates store C4
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fC4, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(0(0.5)1.5)
			//graph export "$results\Figures\5. Marginsplot M1 Violence against state (param).png", replace

			
		//	combomarginsplot fA4 fB4 fC4, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the number of violent events against state forces.", size(small))
			//graph export "$results\Figures\Marginsplot M1 Violence against state.png", replace
			
			coefplot (A4, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (B4, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(C4, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(, glpattern(solid) glwidth(vthin) glcolor(gs15))
			//graph export "$results\Figures\Coefplot M1 Violence against state.png", replace
		
		}
	
		
		* Turf wars
		{
			// Coef plot

			xtreg turf_violence errad_aerea Lturf_violence $controls_weather $controls_economic_performance $controls_distance_airports m_turf_violence i.year if $time_period, fe
			estimates store TU
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fTU, replace)	
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) //ylab(0(20)80)
			// graph export "$results\Figures\6. Marginsplot M1 Turf wars.png", replace
			
		// combomarginsplot fTU, label("Turf wars between armed groups") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the number of violent events between armed groups.", size(small))
			// graph export "$results\Figures\Marginsplot M1 Turf wars.png", replace
			
			coefplot (TU, label(Turf wars between armed groups) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(0(0.00001)0.00003, glpattern(solid) glwidth(vthin) glcolor(gs15))
			// graph export "$results\Figures\Coefplot M1 Turf wars.png", replace
		
		}
		
		
		* Multiple homicide
		{
			// Coefplot
			xtreg percap_total_mulhom errad_aerea Lpercap_total_mulhom $controls_weather $controls_economic_performance $controls_distance_airports mpc_total_mulhom i.year if $time_period, fe
			estimates store A3
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fA3, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(-1(1)6)
			// graph export "$results\Figures\7. Marginsplot M1 Multiple Homicides (overall).png", replace
			
			xtreg percap_insur_mulhom errad_aerea Lpercap_insur_mulhom $controls_weather $controls_economic_performance $controls_distance_airports mpc_insur_mulhom i.year if $time_period, fe
			estimates store B3
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fB3, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(-1(1)6)
			//graph export "$results\Figures\7. Marginsplot M1 Multiple Homicides (insurgents).png", replace

			xtreg percap_param_crimorg_mulhom errad_aerea Lpercap_param_crimorg_mulhom $controls_weather $controls_economic_performance $controls_distance_airports mpc_param_crimorg_mulhom i.year if $time_period, fe
			estimates store C3
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fC3, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(-1(1)6)
			// graph export "$results\Figures\7. Marginsplot M1 Multiple Homicides (param).png", replace
			
		// combomarginsplot fA3 fB3 fC3, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the number of multiple homicides events per hundred thousand" "population.", size(small))
			//graph export "$results\Figures\Marginsplot M1 Multiple Homicides.png", replace
			
			coefplot (A3, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (B3, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(C3, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(0(0.0001)0.0002, glpattern(solid) glwidth(vthin) glcolor(gs15))
			// graph export "$results\Figures\Coefplot M1 Multiple Homicides.png", replace
		
		}

	
/////////////////
		* NEIGHBORS Mechanisms: Displacement of violence
		* Redefining dependent violence variables names
		global mulhom_mechanism_variables "insur_mulhom param_crimorg_mulhom total_mulhom"

		foreach j of global indep_eradication_variables {
			foreach k of global mulhom_mechanism_variables {

				xtreg mpc_`k' `j' Lmpc_`k' $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			//	outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lmpc_`k' $controls_weather $controls_economic_performance $controls_distance_airports) excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes)
			} 
		}
		
		global AoW_mechanisms_variables "total_statev insur_statev param_crimorg_statev turf_violence "

		foreach j of global indep_eradication_variables {
			foreach k of global AoW_mechanisms_variables {

				xtreg m_`k' `j' Lm_`k' $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			//	outreg2 using "$results\Raw tables\M1 `k'.xls", keep(`k' `j' Lm_`k' $controls_weather $controls_economic_performance $controls_distance_airports) excel addtext(Municipalities, All, Municipal FE, Yes, Lagged Y, Yes, Year FE, Yes, Controls, Yes)
			} 
		}
	
	
		* NEIGHBORS: Violence against the state
		{
			// Coef plot
			xtreg m_total_statev errad_aerea Lm_total_statev $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store AM2
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fAM2, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(0(0.5)2)
		//	graph export "$results\Figures\8. Marginsplot M1 Neighbors State Violence (overall).png", replace

			
			xtreg m_insur_statev errad_aerea Lm_insur_statev $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store BM2
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fBM2, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(0(0.5)2)
		//	graph export "$results\Figures\8. Marginsplot M1 Neighbors State Violence (insurgents).png", replace

		
			xtreg m_param_crimorg_statev errad_aerea Lm_param_crimorg_statev $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store CM2
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fCM2, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(0(0.5)2)
		//	graph export "$results\Figures\8. Marginsplot M1 Neighbors State Violence (param).png", replace

			
		//	combomarginsplot fAM2 fBM2 fCM2, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the neighboring municipalities' mean number of violent events" "against state forces.", size(small))
		//	graph export "$results\Figures\Marginsplot M1 Neighbors State Violence.png", replace
			
			coefplot (AM2, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (BM2, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(CM2, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(, glpattern(solid) glwidth(vthin) glcolor(gs15))
	//	graph export "$results\Figures\Coefplot M1 Neighbors State Violence.png", replace
		
		}
		
		
		* NEIGHBORS: Turf wars
		{
			// Coef plot
		
			xtreg m_turf_violence errad_aerea Lm_turf_violence $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store TU1
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fTU1, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) //ylab(0(20)80)
		//	graph export "$results\Figures\9. Marginsplot M1 Neighbors Turf wars.png", replace
			
		//	combomarginsplot fTU1, label("Turf wars between armed groups") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the neighboring municipalities' mean number of violent events" "between armed groups.", size(small))
		//	graph export "$results\Figures\Marginsplot M1 Neighbors Turf wars.png", replace

			coefplot (TU1, label(Turf wars between armed groups) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(-0.00001(0.00001)0.00001, glpattern(solid) glwidth(vthin) glcolor(gs15))
	//	graph export "$results\Figures\Coefplot M1 Neighbors Turf wars.png", replace
	
		}


		* NEIGHBORS: Multiple homicide
		{
			// Coef plot: Neighbors Multiple homicide
			xtreg mpc_total_mulhom errad_aerea Lmpc_total_mulhom $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store AM1
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fAM1, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dknavy)) ylab(-0.5(0.5)1.5)
		//	graph export "$results\Figures\9. Marginsplot M1 Neighbors Multiple Homicides (overall).png", replace
			
			xtreg mpc_insur_mulhom errad_aerea Lmpc_insur_mulhom $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store BM1
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fBM1, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(maroon)) ylab(-0.5(0.5)1.5)
		//	graph export "$results\Figures\9. Marginsplot M1 Neighbors Multiple Homicides (insurgents).png", replace

		
			xtreg mpc_param_crimorg_mulhom errad_aerea Lmpc_param_crimorg_mulhom $controls_weather $controls_economic_performance $controls_distance_airports i.year if $time_period, fe
			estimates store CM1
			margins, at(errad_aerea=(0 30000)) atmeans vsquish saving(fCM1, replace)
			marginsplot, recast(line) recastci(rarea) ciopt(color(gs14)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") /*note("Note: Y axis displays the adjusted linear predictions for the number of violent events per hundred thousand population.", size(small))*/ plotopts(lcolor(dkgreen)) ylab(-0.5(0.5)1.5)
		//	graph export "$results\Figures\9. Marginsplot M1 Neighbors Multiple Homicides (param).png", replace

			
		//	combomarginsplot fAM1 fBM1 fCM1, label("Overall" "Insurgents" "Paramilitaries and Criminal Organizations") recast(line) recastci(rarea) ciopt(color(gs14)) file1opts(color(dknavy)) file2opts(color(maroon)) file3opts(color(dkgreen)) legen(position(6)) xlab(0(10000)30000, format(%9.0fc)) xtitle("Coca crops hectares eradicated via aerial spraying") ytitle("Adjusted linear predictions") title("") note("Note: Y axis displays the adjusted linear predictions for the neighboring municipalities' mean number of multiple" "homicides events per hundred thousand population.", size(small))
		//	graph export "$results\Figures\Marginsplot M1 Neighbors Multiple Homicides.png", replace
			
			coefplot (AM1, label(Overall) msize(medsmall) mcolor(dknavy) ciopts(color(dknavy) lwidth(medthick thick) recast(rcap) msize(large))) (BM1, label(Insurgents) msize(medsmall) mcolor(maroon) ciopts(color(maroon) lwidth(medthick thick) recast(rcap) msize(large)))(CM1, label(Paramilitaries and Criminal Organizations) msize(medsmall) mcolor(dkgreen) ciopts(color(dkgreen) lwidth(medthick thick) recast(rcap) msize(large))), keep(errad_aerea) legend(position(6)) coeflabels(errad_aerea="Aerial Spraying Eradication") vertical yline(0, lcolor(gs6)) ylab(, glpattern(solid) glwidth(vthin) glcolor(gs15))
	//	graph export "$results\Figures\Coefplot M1 Neighbors Multiple Homicides.png", replace
				
		}
}	
